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O ■ ABSTRACT 
(N 
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^2 ! Employing two samples containing of 56 and 59 well-separated FRED (fast 

rise and exponential decay) gamma-ray burst (GRB) pulses whose spectra are 
O I fitted by the Band spectrum and Compton model, respectively, we have inves- 

tigated the evolutionary slope of Ep (where Ep is the peak energy in the uFu 
'^j^'] spectrum) with time during the pulse decay phase. The bursts in the samples 

jnr] I were observed by the Burst and Transient Source Experiment (BATSE) on the 

(-H ' Compton Gamma-Ray Observatory. We first test the Ep evolutionary slope dur- 

0-c_ ing the pulse decay phase predicted by Lu et al. (2007) based on the model of 

O ■ highly symmetric expanding fireballs in which the curvature effect of the expand- 

' ing fireball surface is the key factor concerned. It is found that the evolutionary 

^ ■ slopes are normally distributed for both samples and concentrated around the 

values of 0.73 and 0.76 for Band and Compton model, respectively, which is in 
>■ ' good agreement with the theoretical expectation of Lu et al. (2007). However, 

the inconsistence with their results is that the intrinsic spectra of most of bursts 



in 



■ may bear the Comptonized or thermal synchrotron spectrum, rather than the 
CO ' 

^ ' Band spectrum. The relationships between the evolutionary slope and the spec- 

Q ■ tral parameters are also checked. We show the slope is correlated with Ep of 



O 



time-integrated spectra as well as the photon flux but anticorrelated with the 
lower energy index a. In addition, a correlation between the slope and the in- 
trinsic Ep derived by using the pseudo-redshift is also identified. The mechanisms 
■ of these correlations are unclear currently and the theoretical interpretations are 

^ required. 
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1. INTRODUCTION 

The origin of the gamma-ray burst (GRB) still remains unclear though it has been 
found for more than forty years and many progresses have been made. The best probe of 
early phase of jet attributes and the physical mechanism is the prompt emission, essentially 
the temporal-spectral dependence of GRB pulses due to the absence of some other ways to 
detect them, such as gravitational wave and neutrino detections. The analysis of the GRB 
prompt emission provides us valuable clues to the environment from which the radiation 
emits and the underlying processes giving rise to the phenomenon. The spectral evolution 
is universal and has been studied both over the entire burst, giving the overall behavior and 
over individual pulse structures (see, for instance, the review by Ryde 1999). Pulses are 
common features in a GRB light curve and appear to be the fundamental constituent of it 
(see, e.g., Norris et al. 1996; Stern & Svensson 1996). This single pulse temporal evolution 
is often been described by a fast rise and exponential decay (the so-called FRED shape, see 
Fishman et al. 1994). The spectral hardness decreases during pulse decay and there are two 
empirical relations between the temporal and spectral properties to characterize the spectral 
evolution. One important correlation is that between hardness of the spectrum, the peak 
energy, and fluence (HFC; Liang & Kargatis 1996). The other correlation is that between 
hardness and the instantaneous flux (or intensity) (HIC; Golenetskii et al. 1983). The two 
relations appear to be satisfled by a substantial fraction of GRB pulses during the decay 
phase (see, Crider & Liang 1999; Ryde & Svensson 2000). Most studies concerning these 
correlations examine them in single pulses and do not compare the behavior of pulses within 
a burst. Combining the two correlations, i.e., the HIC and HFC, Ryde & Svensson (2000) 
showed that a power-law HIC and an exponential HFC resulted in the decay phase of the 
pulse following power-law behaviors described by equations 3 and 4 in their paper. Ryde 
& Svensson (2002) studied a sample containing 25 pulses and found that a power law gives 
a better description of the pulse decays than a stretched exponential, the most commonly 
assumed pulse shape so far. They also found that there are no obviously preferred values of 
the power- law index of peak energy, Ep, in the decay phase of pulse. 

Daigne & Mochkovitch (2003) presented a simple, semi-analytical model to interpret 
the GRB temporal and spectral properties in the context of the internal shock model. The 
spectral evolution of synthetic pulses was flrst obtained with standard equipartition assump- 
tions to estimate the post-shock magnetic fleld and the electron Lorentz factor. They found 
that Ep DC t~^, which shows the decay of peak energy during the pulses decay phase follows 
the power law relation with time. They also compared the consistence of the power-law 
index, 6, of their model with the observed index provided by Ryde & Svensson (2002) and 
found that the synchrotron process with standard equipartition assumptions gives a much 
too steep spectral evolution, which showed that another process different from synchrotron 
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one might radiate that energy. 

The observed GRB pulses are beheved to be produced in a relativistically expanding 
and coUimated fireball because of the large energies and the short time-scales involved. The 
so-called Doppler effect (or curvature effect in some paper), which over the whole fireball 
surface would play an important role to account for the observed pulses and spectra of not 
only prompt GRBs but the early X-ray afterglow, is that the photons emitted from the 
regions on the line of sight and off the line of sight are Doppler-boosted by different factors 
and travel different distances before reaching the observer (e.g. Meszaros and Rees 1998; 
Hailey et al. 1999; Qin 2002, Qin et al. 2004, Qin 2008a, 2008b, 2008c). The curvature 
effect model of a cutoff power law spectrum has also been used to model the light curve and 
spectral evolution of the X-ray tail (Zhang et al. 2009). Based on the model of Doppler 
effect Lu et al. (2007) (hereafter Paper I) investigated the evolution of observed spectral 
peak energy Ep and found that the evolutionary curve of Ep undergoes a drop-to- rise-to-decay 
phase. The decay phase of the pulse, where the Doppler effect dominates, always decreases 
monotonically and does not necessarily reflect the corresponding intrinsic spectral evolution. 
They first explored the case the intrinsic spectrum is the Band function (Band et al. 1993) 
for three different local pulses. The three decay phases of Ep within the decay phase of the 
light curve were extracted and performed a linear least-squares fit to the Ep and relative 
observed time, r, and have log Ep — I — 5'logr, with the same slope of 5" = 0.95 ± 0.01 
for the three curves, which indicated that the slope of the decay phase of Ep is not affected 
by the shape of its local pulse. Then they studied the case of intrinsic Comptonized and 
thermal synchrotron spectrum and found that corresponding slope S — 0.75 ± 0.01 for both 
of the two intrinsic spectra, which also showed that the slope is independent of the forms of 
local pulse. 

Although several GRB missions such as HETE-2 and Swift have been launched, the 
Burst and Transient Source Experiment (BATSE; Fishman et al. 1989), aboard the Compton 
Gamma Ray Observatory (CGRO; Gehrels et al. 1994), provided the largest GRB database 
from a single experiment among all the gamma-ray experiments that have detected GRBs. 
The BATSE data are still the most suitable for detailed spectral studies of GRB prompt 
emission, both in quantity and quality. For many of the BATSE GRBs, high time and energy 
resolution data are available. BATSE also provided wider energy coverage than current GRB 
missions. 

In the present work, we select two samples observed by BATSE to investigate the evo- 
lutionary slope of Ep during the decay phase of GRB pulses. Different from what Ryde & 
Svensson (2002) did, we adopt those well-separated FRED pulses compiled by Peng et al. 
(2007). First of all, we check whether the observed evolutionary slope of Ep during the pulse 
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decay phase is consistent with the theoretical predication of Paper I. Examining whether 
the slope of Ep during the decay phase of pulse is related to other parameters such as Ep 
of time-integrated spectra is another motivation of this work. In section 2, we present the 
sample description and spectral modeling. The results are given in section 3. Discussion and 
conclusions are presented in the last section. 

2. SAMPLE DESCRIPTION AND SPECTRAL MODELING 

In this work we only select the GRB sample compiled by Peng et al. (2007) consisting 
of two samples provided by Kocevski et al. (2003) and Norris et al. (1999), in which the 
bursts are found to contain individual FRED pulses. The data with durations longer than 
2 s are provided by the BATSE's LADs instruments on board the CGRO spacecraft, which 
provide discriminator rate with 64 ms resolution from 2.048 s before the burst to several 
minutes after the trigger (Fishman et al. 1994) (for more details of the sample selection, see 
Kocevski et al. 2003. and Norris et al. 1999). In order to obtain the evolutionary slope 
of Ep during the pulse decay phase (we use the symbol S to denote the evolutionary slope 
of Ep during the pulse decay phase) we must study the time-resolved spectroscopy. The 
time-resolved spectral analysis requires data with a high enough signal-to-noise (S/N) ratio 
to obtain better statistical results. Similar to Ryde & Svensson (2002) and Firmani et al. 
(2008) wc also choose those bursts whose peak flux, Fp, above 1.8 photons cm^^s^^ on a 256 
ms timescale to ensure a reasonable determination of the spectral parameters. In this way 
78 pulses meet the criteria. 

Before performing the spectral modeling we must select the type of detector, data type, 
and time and energy interval. Following the Kaneko et al. (2006) (hereafter Paper II) we also 
only use the LAD data mainly to take advantage of its larger effective area. Three LAD data 
types used in this work are Medium Energy Resolution data (MER), High Energy Resolution 
Burst data (HERB), and Continuous data (CONT) (for more detailed description about the 
characteristics of each data type one can refer to Paper II) . To perform detailed time- resolved 
spectroscopy it has been shown by Preece et al. (1998) that a S/N ~45 is needed. However, 
a high S/N will lead to time- resolved spectra consisting of only a few broader time bins. 
Our present work focus on the analysis of the S. Therefore the aim of our spectral analysis, 
of every time bin, is mainly to determine the Ep and allowing a deconvolution of the count 
spectrum to flnd the energy spectrum. In order to arrive at reliable results wc need as many 
time bins as possible to study the S. For this purpose, we adopt the criterion that the decay 
phase of the pulses should have at least five time bins with 5"/ A?" > 30 to be included in the 
study. We check that the results are consistent with higher S/N ratios. This gives us the 
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possibility to study the burst pulses with higher time resolution, especially for the later time 
bins, which is of great importance for our study. 

Similar to Paper II we also use the spectral analysis software RMFIT, which was specif- 
ically developed for burst data analysis by the BATSE team (Mallozzi et al. 2005). It in- 
corporates a fitting algorithm MFIT that employs the forwardfolding method (Briggs 1996), 
and the goodness of fit is determined by minimization. One advantage of MFIT is that it 
utilizes model variances instead of data variances, which enables more accurate fitting even 
for low-count data (Ford et al. 1995). 

We analyze both time-resolved and time-integrated spectra for each pulse of our sample. 

In addition, the entire pulse must be included in our analysis since these pulses are well- 
separated. The time-resolved and time-integrated spectra are modeled with two photon 
models. One is the most used so-called Band function (Band et al. 1993): 

, _^f(li)"«P(-^) ^<^^c, 

" I (^^S^)"-''exp(/J - .)(-^Y E > B„ 

where, Ec = (a — P) ^"^^ = (a — /9)£'o- The model consists of four parameters: the amplitude 
A in photons s~^ cm~^ keV~^, a low-energy spectral index «, a high energy spectral index 
/3, and a vFy peak energy -EpeaA- in kcV, which is related to the c-folding energy, Eq. 

Paper II have showed that the Comptonizcd Model (COMP model) tends to be prefer- 
able in fitting time-resolved spectra due to better values compared with that of BAND 
model. Additionally, many BATSE GRB spectra lack high-energy photons (Pendleton et al. 
1997). Those no-high-energy spectra are usually fitted well with this model. This also shows 
the existence of a number of spectra without high-energy component. Therefore, the spectra 
are also modeled with the COMP model. It is a low-energy power law with an exponential 
high-energy cutoff, which is equivalent to the BAND model without a high-energy power 
law, namely /3 — >cx:, and has the form 

JcoMpiE) = Ai^reM-^^). (2) 

Epiy was always fixed at 100 keV , therefore, the model consists of three parameters: A, a, 
and Epeak- 

We always chose the data taken with the detector that was closest to the line of sight to 
the location of the GRB, as it has the strongest signal. For the case of weaker bursts the MER 
data is preferable because it has much finer time resolution. Whereas for the bright bursts 
the HERB data is preferable due to its higher energy resolution. A background estimate 
is made using the High Energy Resolution data (HER) data for the HERB data, covering 
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several thousand seconds before the trigger and after the HERB accumulation is finished. In 
the MER data, the CONT data are used as background. The hght curve of the background, 
during the outburst, is modeled with a second- or third-order polynomial fit. The usable 
energy range for spectral analysis is ~ 30 keV - ~ 2 MeV for all the bursts. The lowest seven 
channels of HERB and two channels of MER and CONT are usually below the electronic 
lower energy cutoff and are excluded. The highest few channels of HERB and normally the 
very highest channel of MER and CONT are unbounded energy overflow channels and also 
not usable (see. Paper H). The procedures give us all the spectral parameters modeled by 
BAND and COMP model for every time bin of the whole pulses besides energy flux. Similar 
to Peng et al. (2009) the following data points are excluded: 

(1) the resulting per degree of freedom is -1 when fitting each time-resolved spectrum 
because in this case the nonlinear fitting of corresponding time-resolved spectra is failure. 

(2) q; < -2 and /? > -2 for the BAND model, a < -2 for the COMP model. As Paper 
n pointed out that, the fitted Ep^ak represents the actual peak energy of the vFy spectrum 
only in the case of a > -2 and /3 < -2 for the BAND model and a > -2 for the COMP 
model. 

(3) the uncertainty of its corresponding Ep^ak is larger 50% than itself. In this way, we 
can obtain the best statistic. 

With these criterions, we can obtain useful data points of -Epeafe as well as the other 
parameters. At the end, we have 56 and 59 pulses for Band model and Comp model, which 
are denoted with sample 1 and 2, respectively. 

We focus, in this study, our attentions on the decay phase of pulses. Hence the initial 
time of the decay phase must be determined. Similar to Peng et al. (2006, 2009) the 
subtract-background pulses, combining the data from the BATSE four channels are modeled 
with the function presented in equation (22) of Kocevski et al. (2003) because we find that 
this function can well describe the observed profile of a FRED pulse. In addition, a fifth 
parameter, t^, which measures the offset between the start of the pulse and the trigger time, 
is introduced. 

where tm is the time of the pulse's maximum flux, F^; r and d are the power-law rise and 
decay indexes, respectively. Note that equation (1) holds for t > —to, when t < —to we take 
F(t) = 0. 

To obtain an intuitive view of the result of the fit, we develop and apply an interactive 
IDL routine for fitting pulses in bursts, which allows the user to set and adjust the initial 
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pulse parameters manually before allowing the fitting routine to converge on the best-fitting 
model via the reduced minimization. The fits are examined many times to ensure that 
they arc indeed the best ones. The tm of each pulse are found and served as the starting 
time of the pulse decay phase. The values, tm, io fitting chi-square are listed in Tables 
1 and 2. In this way the data points of Ep within the decay portions of the FRED pulses 
are extracted. 



3. analysis result 

3.1. the evolutionary slope seen in the FRED pluses 

Since Paper I have shown that the evolutionary slope in pulse decay phase can be well 
modeled with a single power-law form we first adopt this form to fit the temporal evolution 
of Ep during the decay phase of pulse. 

m - JJ^, (4) 



The example plots of the fits within the decay phases are given in Figures 1 and 2, which 
correspond to samples 1 and 2, respectively. 

In order to compare these results with the theoretical ones presented by Paper I, we also 
plot them in Figures 1 and 2 in a logarithmic manner for samples 1 and 2, respectively. In 
addition, the evolutionary slope of the fitted peak energy Ep should be dlogEp/ d\og{t -\-tQ) 
rather than dlogEp/ d\og{t -\- ttrigger), where to is the starting point of the rising segment 
of the FRED pulse. For each pulse of our sample we use LINFIT (Press 1992) to perform 
a linear least-squares fit at a 1 a confidence level to the two quantities and have log Ep = 
W — S'log(t + to)- The evolutionary slope, S, are listed in Tables 1 and 2. We first check the 
consistence between the power-law index 6 and the slope S. Figure 3 (right panel) gives the 
scatter plot of them, which shows they are indeed consistent. 

The distribution histogram of S is indicated in Figure 3. We find from Figure 3 for 
the two samples that: (i) the slopes vary from 0.23 to 1.66 for sample 1 and from 0.35 to 
1.47 for sample 2; (ii) the distributions of slopes arc normal and the corresponding mean 
values and standard deviations are around 0.73 and 0.22 for sample 1 and 0.76 and 0.22 for 
sample 2, respectively; (iii) the corresponding values of the medians are 0.74 and 0.76 for 
samples 1 and 2, respectively. These indicate that for most of the pulses the slopes are good 
consistent with that predicted by paper I, which showed that the evolutionary slopes of the 
pulse decay phase are 5"! = 0.95 ± 0.01 and S2 — 0.75 ± 0.01 corresponding to the intrinsic 
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Fig. 1. — The example plots of spectral and temporal behaviors of the pulses in sample 1. 
The top panels show the evolutions of the peak energy Ep of the time-resolved spectra (left 
axes) and their corresponding pulse light curves (all four energy channels are used) (right 
axes). The best fit of pulse is indicated with a long dashed curve, the vertical dotted hue 
represents the peak time of pulses, while the dotted curve is the fitting line of the Ep using 
a single power-law function. The bottom panel indicates the developments of the Ep within 
the decay phase of the pulses and the long dashed lines represent the best linear fitting. 
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Fig. 2. — The example plots of spectral and temporal behaviors of the GRB pulses in sample 
2. The meaning of the symbols are the same as figure 1. 
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Band function spectrum and Comptonized or thermal synchrotron spectrum, respectively. 
Therefore, within the limit of uncertainty the corresponding intrinsic spectra of most of 
bursts (about 25 for sample 1 and 29 for sample 2) may bear the intrinsic Comptonized or 
thermal synchrotron spectrum, whereas only a small fraction of bursts (about 16 for sample 
1 and 12 for sample 2) may be associated with the intrinsic Band spectra. There are also 
some pulses (about 15 for sample 1 and 18 for sample 2) whose slopes are inconsistent with 
both theoretical values, which is similar to the result investigated by Paper I. 

3.2. the relation between evolutionary slope and the spectral parameters 

We suspect there is a relation between the first Ep (we denote it as Ep ^ax) of the decay 
pulse and the slope because wc find that: (i) the first Ep is usually the maximal value of 
the decay pulse; (ii) the larger the first Ep is, the steeper of the fitting curves is. Figure 4 
demonstrates the relation between Ep^^nax and S for the the samples 1 (panel (a)) and 2 (panel 
(b)). An correlation between the two quantities is identified for both samples. For the sample 
1 a finer regression analysis (log 5" — A-\- B x logEp) yields the intercept A — —1.78 ± 0.08 
and the slope B = 0.73±0.03, with the correlation coefficient R = —0.67 {N = 56,p < 10~^). 
While for the sample 2, the analysis produces the intercept A = —1.66 ± 0.06 and the slope 
B = 0.66 ± 0.02, with the correlation coefficient R ^ -0.68 {N = 59,p < 10"^). 

Sakamoto et al. (2008) have showed that the Ep is correlated with the photon index 
(r) derived from a simple power-law model. Motivated by this, we also examine if there 
are also correlation between the Ep and the S since the S corresponds to the power-law 
index 6. Figure 5 displays the Ep of time-integrated spectra versus S. A clear correlated 
relationship between the two quantities is found. A straight line is fitted to the points: 
(i) logS = -1.16(±0.032) - 0.50(±0.01)logEp, with the correlation coefficient R = 0.65 
(N=56,p < 10-^) for sample 1; (ii) S = -1.22(±0.022) -0.51(±0.01)Ep, with the correlation 
coefficient R = 0.67 (N= 59, p < 10"^) for sample 2. 

From the Figures 4 and 5 we find that both the Ep^^nax and the Ep of time-integrated 
spectra are correlated with the S. We deduce that there must be some relation between 
them. Figure 6 gives the scatter plot of the Ep^rnax versus the Ep of time-integrated spectra, 
the corresponding correlation coefficient it! = 0.94 (N= 56, p < 10~^) for the sample 1 and 
R — 0.92 (N= 59, p < 10~^) for sample 2, which indicates there are indeed strong correlation 
between them. 

We also examine the relation between the S and lower energy index a as well as photon 
flux, /, (the energy range is 25-1800 keV). Figure 7 shows the plot of the a versus S. A 
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Fig. 3. — The distribution histograms of slope, S, (left panel) for sample 1 (solid line) and 
sample 2 (dashed line), where the curves represent the Gaussian fit to the two distributions. 
The right panel is the scatter plot of S vs. power- law index 6 for sample 1. 
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Fig. 4. — The evolutionary slope S vs. the Ep^max of decay pulse for sample 1 (a) and sample 
2 (b), where the long dashed lines represent the best fitting lines. 
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Fig. 5. — The log S vs. logEp for sample 1 (a) and sample 2 (b). 




Fig. 6. — Ep of timc-intcgratcd spectra vs. the Ep ^ax ior sample 1 (a) and sample 2 (b). 
The strong correlation is indicated for both samples. 
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anticorrelation is suggested for both of samples with the correlation coefficient R — —0.57 
(N= 56, p < 10-^) for the sample 1 and i? = -0.44 (N= 59, p = 4.15 x 10"^) for sample 2. 

Figure 8 plots the / versus S. It is found that S is also correlated with /. The regression 
analysis give the best fitting line: (i) log ,5 = -0.26(±0.01) + 0.27(±0.01) log/, with the 
correlation coefficient R = 0.82 (N=56, p < 10"^) for sample 1; (ii) log-S = -0.22(±0.02) + 
0.24(±0.01) log/, with the correlation coefficient R = 0.80 (N= 59, p < 10"^) for sample 2. 

Since the S is related to the Ep, we wonder if there are also correlation between the S 
and intrinsic peak energy £'p j. To verify the relation we need to estimate the redshifts of our 
selected bursts due to unknown redshift information of BATSE bursts. We attempt to use 
the pseudo-redshift derived by Yonetoku et al. (2004). This pseudo-redshifts were estimated 
based on a new relation between the spectral peak energy Ep and the 1 s peak luminosity, 
which was derived by combining the data of Ep and the peak luminosities by BeppoSAX 
and BATSE. In addition, it looks considerably tighter and more reliable than the relations 
suggested by the previous works (e.g. Amati et al. 2002). The intrinsic peak energy, Ep^i, 
in the rest frame is related with the observed Ep in the observer frame by Ep^i = Ep{l + z) . 
There are 49 and 53 pulses with pseudo-redshift information in samples 1 and 2, respectively. 
The estimated relation between them is displayed in Figure 9. 

An correlated relationship between the two quantities is also identified for both samples. 
The best fit functional form of this relation is log S = -0.923(±0.047)-0.298(±0.017) log^;^,^, 
with the correlation coefficient R — 0.67 (N=50, p < 10~^) for sample 1. For the sample 2 
logS = -1.388(±0.061) - 0.492(±0.022)log£;p,i, with the correlation coefficient R = 0.50 
(N= 53, p < 10~^). Comparing it with the case of observed relation between S and Ep the 
intrinsic relation seems more loosely defined for sample 2. 
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Fig. 7. — Evolutionary slope S vs. low energy index a for sample 1 (a) and sample 2 (b), 
where the long dashed lines represent the best fitting lines. 
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Fig. 8. — Evolutionary slope S vs. photon flux / for sample 1 (a) and sample 2 (b), where 
the long dashed lines represent the best fitting lines. 



Table 1. A list of burst sample and various parameters for sample 1. 
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0.72 


it 


0.21 


0.25 


999 


11.55 ± 0.31 


3.99-6.00 


9 


-3.22 ± 0.14 


3.99 ± 0.00 


2.37 


381.13 it 33.72 


-0.94 it 0.05 


1.05 


± 


0.33 


1.18 


it 


0.41 


0.51 


1406 


1.97 ± 0.13 


3.35-20.00 


10 


0.95 ± 0.07 


3.35 ± 0.05 


1.07 


118.71 it 6.31 


-0.91 it 0.11 


0.65 


± 


0.13 


0.62 


it 


0.13 


0.63 


1883 


5.20 ± 0.18 


1.31-7.00 


6 


0.25 ± 0.07 


1.31 ± 0.02 


0.97 


249.72 it 19.61 


-1.372 ± 0.05 


0.94 


± 


0.27 


1.01 


it 


0.20 


0.56 


1956 


2.57 ± 0.13 


2.91-15.00 


8 


0.07 ± 0.29 


2.91 ± 0.04 


1.21 


126.23 it 8.82 


-1.17 it 0.11 


1.11 


± 


0.14 


1.00 


it 


0.17 


0.95 


1989 


2.73 ± 0.15 


116.54-130.00 


5 


-111.58 ± 0.20 


116.54 ± 0.05 


1.32 


89.32 it 4.93 


-0.99 ± 0.15 


0.80 


± 


0.05 


0.84 


it 


0.09 


1.11 


2083- 1 


45.42 ± 0.46 


1.05-6.00 


17 


0.36 ± 0.22 


1.04 ± 0.00 


5.98 


363.25 ± 8.24 


-1.40 ± 0.02 


1.24 


± 


0.09 


1.42 


± 


0.03 


1.15 


2083- 2 


45.42 ±0.46 


8.68-20.00 


42 


-4.02 ± 0.78 


8.68 ± 0.01 


2.26 


69.95 it 1.40 


-1.31 ± 0.07 


0.83 


it 


0.34 


0.84 


it 


0.04 


0.79 


2138 


7.00 ± 0.20 


1.29-15.00 


11 


1.28 ± 0.21 


1.29 ± 0.06 


1.17 


205.92 it 17.4 


-0.49 it 0.18 


0.59 


it 


0.02 


0.49 


it 


0.15 


0.5801 


2193 


1.55 ± 0.13 


10.32-42.00 


15 


-0.04 ± 0.5 


10.32 ± 0.24 


0.98 


274.11 it 6.48 


0.38 it 0.06 


0.74 


it 


0.46 


0.78 


it 


0.05 


1.08 1 


2387 


3.86 ± 0.16 


6.49-40.00 


38 


-0.53 ± 0.27 


6.49 ± 0.07 


1.15 


132.22 ± 4.17 


-0.27 ± 0.09 


0.76 


± 


0.04 


0.66 


± 


0.04 


1.48 


2484 


1.55 ± 0.13 


2.01-10.00 


6 


2.66 ± 0.85 


2.01 ± 0.07 


1.09 


107.55 ± 1.60 


-0.70 ± 0.06 


0.67 


± 


0.13 


0.86 


± 


0.15 


0.85 


2519 


1.53 ± 0.14 


0.59-15.00 


13 


0.45 ± 0.02 


0.59 ± 0.06 


1.02 


163.76 ± 20.23 


-0.77 ± 0.22 


0.67 


± 


0.07 


0.60 


± 


0.07 


1.00 


2662 


1.52 ± 0.14 


1.31-12.00 


6 


0.30 it 0.05 


1.31 ± 0.07 


1.11 


149.64 it 11.03 


-0.90 it 0.15 


0.53 


it 


0.07 


0.62 


it 


0.09 


1.15 


2665 


1.99 ± 0.15 


1.42-8.00 


7 


1.22 ± 0.32 


1.42 ± 0.05 


0.89 


105.03 it 2.25 


-0.15 it 0.13 


0.73 


± 


0.22 


0.74 


it 


0.25 


0.50 


2700 


4.06 ± 0.18 


53.7-61.00 


14 


-52.84 ± 0.01 


53.71 ± 0.03 


1.37 


242.36 it 15.13 


-1.33 it 0.04 


0.61 


it 


0.21 


0.54 


± 


0.14 


0.12 


2880 


2.90 ± 0.14 


0.47-3.00 


6 


0.34 ± 0.04 


0.47 ± 0.01 


1.43 


124.33 it 6.46 


-0.74 it 0.16 


0.46 


it 


0.12 


0.50 


it 


0.14 


0.40 


2919 


5.77 ± 0.19 


0.34-7.00 


9 


1.57 ± 0.12 


0.34 it 0.02 


1.39 


380.43 it 36.64 


-1.28 ± 0.04 


0.79 


± 


0.19 


0.86 


± 


0.29 


0.82 


3003 


2.83 ± 0.16 


9.68-24.00 


10 


8.94 ± 7.15 


9.68 it 0.08 


1.02 


472.83 it 49.82 


-1.22 ± 0.04 


0.93 


± 


0.22 


1.07 


± 


0.73 


0.40 


3143 


2.59 ± 0.14 


0.69-4.00 


8 


0.64 ± 0.25 


0.69 it 0.02 


0.96 


99.09 it 16.27 


-0.85 ± 0.39 


0.77 


± 


0.44 


0.80 


± 


0.39 


1.43 


3256 


1.76 ± 0.11 


1.38-9.00 


20 


0.55 ± 0.08 


1.38 it 0.07 


1.06 


142.93 it 7.39 


-0.50 it 0.15 


0.65 


it 


0.05 


0.70 


it 


0.07 


0.90 


3257 


3.06 ± 0.13 


3.51-35.00 


27 


0.60 ± 0.06 


3.51 it 0.06 


1.10 


196.33 it 5.15 


-0.24 it 0.06 


0.45 


it 


0.02 


0.57 


it 


0.04 


1.50 


3290 


10.71 ± 0.18 


2.98-4.00 


5 


-2.73 ± 0.00 


2.98 it 0.01 


2.52 


50.84 it 10.53 


-1.69 it 0.21 


0.41 


± 


0.30 


0.45 


it 


0.16 


0.60 


3415 


9.16 ± 0.19 


11.56-15.00 


5 


-10.84 ± 0.02 


11.56 it 0.01 


1.44 


245.74 it 28.54 


-1.21 it 0.07 


0.61 


it 


0.33 


0.70 


it 


0.30 


1.83 


3648-2 


5.70 ± 0.15 


23.85-36.00 


5 


-19.47 ± 0.19 


23.85 it 0.18 


0.97 


75.22 ± 3.68 


0.73 ± 0.45 


0.52 


it 


0.08 


0.55 


± 


0.10 


1.50 


3648-3 


5.70 ± 0.15 


41.04-48.00 


12 


-32.20 ± 1.44 


41.04 ± 0.02 


1.36 


173.94 ± 6.30 


-0.81 ± 0.05 


1.49 


± 


0.02 


1.62 


± 


0.09 


0.85 


3765 


25.29 ± 0.27 


66.15-73.00 


12 


-63.70 ± 0.24 


66.15 ± 0.01 


1.45 


294.73 ± 7.20 


-1.09 ± 0.02 


1.03 


± 


0.19 


1.07 


± 


0.08 


2.05 


3870 


13.90 ± 0.23 


0.50-6.00 


9 


0.46 ± 0.38 


0.50 it 0.01 


1.32 


118.03 it 19.34 


-1.61 it 0.09 


0.76 


it 


0.06 


0.72 


it 


0.14 


0.72 
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Fig. 9. — The evolutionary slope S vs. intrinsic peak energy Ep^i for sample 1 (a) and sample 
2 (b), where the long dashed lines represent the best fitting lines. 



Table 1 — Continued 



Trigger 


(phs cm 


s ) 


• 

time interval 
(V 


n 


(V 


(.sj 


x; 




TP 

\KC V ) 

















Q 






X 


29 


3875 


2.80 


± 


0.11 


0.27-3.00 


8 


0.14 ± 0.01 


0.27 ± 0.02 


1, 


.61 


69.64 ± 4.02 


-1.23 


± 


0.25 


0.47 


± 


0.21 


0.36 


± 


0. 


19 


1. 


,14 


3954 


8.19 


± 


0.19 


0.78-6.00 


8 


2.26 ± 1.47 


0.78 ± 0.01 


1, 


.05 


297.13 ±31.34 


-1.14 


± 


0.06 


1.05 


± 


0.08 


0.98 


± 


0. 


27 


0, 


,86 


4350 


3.27 


± 


0.12 


14.05-20.00 


8 


-13.08 ± 0.03 


14.05 ± 0.04 


1, 


.49 


130.96 ± 11.84 


-1.55 


± 


0.08 


0.56 


± 


0.14 


0.66 


± 


0, 


,17 


1, 


10 


5478 


2.96 


± 


0.12 


2.08-12.00 


11 


0.52 ± 0.05 


2.08 ± 0.04 


0, 


.97 


156.36 ± 4.63 


-0.02 


± 


0.09 


0.70 


± 


0.12 


0.75 


± 


0, 


,08 


1, 


,58 


5517 


1.77 


± 


0.11 


0.83-6.00 


9 


0.18 ± 0.44 


0.83 ± 0.07 





.97 


157.15 ± 13.86 


-1.18 


± 


0.13 


0.73 


± 


0.21 


0.77 


± 


0. 


,21 


0. 


,62 


5523 


3.67 


± 


0.14 


1.05-5.00 


5 


0.33 ± 0.18 


1.05 ± 0.03 


1, 


.16 


171.85 ± 11.44 


-1.01 


± 


0.09 


1.01 


± 


0.22 


1.10 


± 


0. 


,25 


0. 


,67 


5601 


4.49 




0.14 


7.67-15.00 


25 


-5.04 ± 0.18 


7.67 ± 0.02 


1, 


.31 


205.55 ±9.44 


-0.72 


± 


0.06 


1.06 


± 


0.05 


1.02 


± 


0. 


,11 


0. 


,99 


6159 


1.94 


± 


0.12 


3.17-12.00 


7 


2.28 ± 8.25 


3.17 ± 0.11 


1 


.01 


68.91 ± 5.12 


-0.09 


± 


0.56 


0.65 


± 


0.17 


0.70 


± 


0. 


,21 


0. 


,89 


6397 


5.78 




0.15 


3.42-25.00 


13 


0.55 ± 0.05 


3.42 ± 0.02 


1, 


.40 


192.55 ± 5.44 


-0.65 


± 


0.04 


0.55 


± 


0.06 


0.64 


± 


0, 


,06 


1, 


,72 


6504 


2.33 


± 


0.12 


3.09-17.00 


10 


1.04 ± 0.51 


3.09 ± 0.07 


0, 


.98 


155.32 ± 4.29 


0.43 


± 


0.11 


0.74 


± 


0.08 


0.77 


± 


0, 


,06 


0, 


,62 


6621 


6.71 


± 


0.16 


32.53-37.00 


10 


-29.56 ± 0.48 


32.53 ± 0.02 


1, 


.09 


95.04 ± 4.87 


-0.98 


± 


0.13 


0.89 


± 


0.23 


0.81 


± 


0. 


13 


2, 


,58 


6625 


1.81 


± 


0.13 


5.22-23.00 


13 


1.51 ± 0.84 


5.22 ± 0.11 


1, 


.07 


82.74 ± 1.98 


-0.89 


± 


0.11 


0.79 


± 


0.10 


0.76 


± 


0. 


,11 


0. 


,86 


6657 


1.86 


± 


0.21 


4.26-47.00 


5 


1.31 ± 0.11 


4.26 ± 0.19 


1, 


.06 


70.41 ± 11.14 


1.03 


± 


1.38 


0.36 


± 


0.12 


0.41 


± 


0. 


09 


0. 


,17 


6930 


5.54 


± 


0.18 


31.83-35.00 


9 


-29.94 ± 0.15 


31.83 ± 0.02 


1 


.06 


95.18 ± 4.34 


-0.64 


± 


0.13 


1.12 


± 


0.17 


1.02 


± 


0. 


16 


0, 


19 


7293 


2.95 


±0.11 


3.66-40.00 


19 


0.46 ± 0.05 


3.66 ± 0.06 


1, 


.27 


161.44 ± 3.26 


0.32 


± 


0.08 


0.60 


± 


0.02 


0.50 


± 


0, 


,02 


0, 


,99 


7295 


3.26 


± 0.17 


2.23-9.00 


6 


1.16 ± 0.64 


2.23 ± 0.06 


1, 


.02 


344.84 ± 14.25 


0.49 


± 


0.10 


0.53 


± 


0.12 


0.61 


± 


0, 


13 


0, 


,72 


7475 


3.69 


± 0.14 


9.15-30.00 


16 


-0.05 ± 0.11 


9.15 ± 0.05 


1, 


.84 


173.14 ± 12.34 


-1.29 


± 


0.05 


0.64 


± 


0.04 


0.69 


± 


0. 


,09 


0. 


,67 


7548 


2.95 


± 0.12 


3.77-6.00 


5 


-1.64 ± 0.22 


3.77 ± 0.03 


0, 


.98 


217.14 ± 14.34 


-0.91 


± 


0.07 


1.14 


± 


0.71 


1.05 


± 


0. 


,78 


0. 


,15 


7588 


2.08 


± 0.12 


2.74^13.00 


5 


0.44 ± 0.27 


2.74 ± 0.07 


0, 


.84 


78.84 ± 2.35 


-0.68 


± 


0.17 


0.44 


± 


0.08 


0.54 


± 


0. 


,10 


0. 


,34 


7638 


1.75 


± 0.11 


1.10-9.00 


16 


0.92 ± 0.63 


1.10 ± 0.05 


1, 


.47 


61.34 ± 3.23 


-0.92 


± 


0.51 


0.32 


± 


0.11 


0.41 


± 


0. 


,11 


0. 


,34 


7648 


1.55 


± 0.10 


4.43-31.00 


10 


0.89 ± 0.46 


4.43 ± 0.14 


1, 


.02 


252.93 ± 9.04 


-0.44 


± 


0.06 


0.58 


± 


0.08 


0.57 


± 


0, 


,08 


0, 


,40 


7711 


3.66 


± 0.13 


1.98-8.00 


10 


2.01 ± 1.78 


1.98 ± 0.04 


1, 


.28 


211.03 ± 12.83 


-1.24 


± 


0.05 


1.06 


± 


0.19 


1.05 


± 


0, 


18 


0, 


,98 


8049- 1 


1.59 


± 0.09 


8.03-20.00 


7 


15.45 ± 9.02 


8.03 ± 0.16 


1, 


.02 


149.13 ± 5.45 


-0.32 


± 


0.11 


1.21 


± 


0.55 


1.02 


± 


0. 


,57 


0, 


,37 


8049-2 


1.59 


± 0.09 


30.99-80.00 


9 


13.15 ± 10.17 


30.99 ± 0.18 


0, 


.96 


115.03 ± 4.38 


-1.26 


± 


0.08 


1.17 


± 


0.15 


1.04 


± 


0. 


,17 


1. 


,68 



Note. — Fp denotes the peak flux on a 256 ms timescale; time interval is the interval of start and end time of decay phase of selected pulses. ^iid are 
the fitting xt of pulse light curves and the power- law decay of Ep, respectively. 



Table 2. A list of burst sample and various parameters for sample 2. 



Trigger Fp time interval 


n 


to 




Ep 


a 


S 


S 


X?2 


(phs cm~'^s~^) (s) 




(s) 


(s) 


(keV) 










Band photon model 



563 


1.89 ± 0.14 


1.84-32.00 


12 


0.66 ± 0.04 


1.84 ± 0.05 


1, 


.25 


169.53 ± 4.74 


-0.58 


it 


0.07 


0.75 


it 


0.04 


0.78 


it 


0.09 


1, 


,62 


907 


3.57 ± 0.17 


1.67-16.00 


16 


0.31 ± 0.04 


1.67 ± 0.03 


1, 


.38 


198.34 ± 3.30 


0.26 


it 


0.06 


0.67 


it 


0.02 


0.74 


it 


0.01 


1, 


,61 


914 


2.53 ± 0.16 


0.61-5.00 


6 


1.17 ± 0.45 


0.61± 0.02 


0. 


.92 


116.46 ± 4.77 


-1.29 


it 


0.08 


1.20 


it 


0.21 


1.21 


it 


0.20 


1, 


,25 


973_1 


5.29 ± 0.20 


2.79-22.00 


21 


-0.41 ± 0.07 


2.79 ± 0.03 


1, 


.34 


412.52 ± 14.11 


-1.26 


it 


0.02 


0.95 


it 


0.09 


1.00 


it 


0.06 


1, 


,41 


973-2 


5.29 ± 0.20 


24.09-35.00 


7 


-22.14 ± 0.01 


24.09 ± 0.03 


1, 


,24 


296.21 ± 24.92 


-1.51 


± 


0.04 


0.82 


it 


0.15 


0.78 


it 


0.14 


1, 


,25 


999 


11.55 ± 0.31 


3.99-6.00 


13 


-3.22 ± 0.14 


3.99 ± 0.00 


2, 


,37 


600.93 ± 31.72 


-1.12 


± 


0.03 


1.31 


it 


0.31 


1.34 


it 


0.26 


0, 


,98 


1406 


1.97 ± 0.13 


3.35-37.00 


14 


0.95 ± 0.07 


3.35 ± 0.05 


1, 


,07 


153.71 ± 4.84 


-1.26 


± 


0.05 


0.79 


it 


0.09 


0.81 


it 


0.08 


0, 


,78 


1467 


2.26 ± 0.13 


4.47-28.00 


5 


0.69 ± 0.52 


4.47 ± 0.04 


0, 


,99 


111.64 ± 3.73 


-1.01 


± 


0.12 


0.41 


it 


0.04 


0.44 


it 


0.04 


1, 


,81 


1733 


3.00 ± 0.15 


3.46-30.00 


5 


6.26 ± 3.25 


3.46 ± 0.04 


1, 


,41 


649.54 ± 79.23 


-1.17 


± 


0.04 


1.56 


it 


0.34 


1.21 


it 


0.17 


1, 


,66 


1883 


5.20 ± 0.18 


1.31-12.00 


7 


0.25 ± 0.07 


1.31 ± 0.02 


0, 


,97 


248.32 ± 16.72 


-1.37 




0.04 


0.68 


± 


0.17 


0.74 


it 


0.15 


1, 


,56 


1956 


2.57 ± 0.13 


2.91-15.00 


8 


0.07 ± 0.29 


2.91 ± 0.04 


1, 


,21 


146.23 ± 6.24 


-1.33 




0.06 


0.76 


it 


0.14 


0.69 


it 


0.13 


1, 


,65 


1989 


2.73 ± 0.15 


116.54-130.00 


5 


-111.58 ± 0.20 


116.54 ± 0.05 


1, 


,32 


105.32 ± 3.29 


-1.33 


± 


0.06 


0.93 


it 


0.05 


1.00 


it 


0.09 


1, 


,01 


2083- 1 


45.42 ± 0.46 


1.05-6.00 


17 


0.36 ± 0.22 


1.04 ± 0.00 


5, 


,98 


404.25 ± 6.28 


-1.06 


it 


0.01 


1.34 


it 


0.09 


1.41 


it 


0.03 


1, 


,15 


2083- 2 


45.42 ± 0.46 


8.68-20.00 


42 


-4.02 ± 0.78 


8.68 ± 0.01 


2. 


,26 


64.31 ± 2.54 


-1.64 


it 


0.03 


0.85 


it 


0.34 


0.78 


it 


0.04 


0, 


,77 


2138 


7.00 ± 0.20 


1.29-15.00 


17 


1.28 ± 0.21 


1.29 ± 0.06 


1. 


,17 


205.92 ± 17.4 


-0.49 


it 


0.18 


0.45 


it 


0.02 


0.58 


it 


0.15 


0, 


,58 


2193 


1.55 ± 0.13 


10.32-42.00 


15 


-0.04 ± 0.5 


10.32 ± 0.24 


0, 


,98 


295.31 ± 15.18 


0.23 


it 


0.07 


0.75 


it 


0.06 


0.96 


it 


0.04 


1, 


,25 


2387 


3.86 ± 0.16 


6.49-40.00 


15 


-0.53 ± 0.27 


6.49 ± 0.07 


1, 


,15 


167.22 ± 3.13 


-0.75 


it 


0.04 


0.75 


it 


0.03 


0.77 


it 


0.03 


2, 


,48 


2484 


1.55 ± 0.13 


2.01-30.00 


8 


2.66 ± 0.85 


2.01 ± 0.07 


1, 


,09 


106.55 ± 1.65 


-0.72 


it 


0.03 


1.05 


it 


0.13 


0.97 


it 


0.11 


0, 


,99 


2519 


1.53 ± 0.14 


0.59-15.00 


7 


0.45 it 0.02 


0.59 ± 0.06 


1, 


,02 


192.76 ± 14.73 


-0.97 


± 


0.12 


0.75 


it 


0.07 


0.75 


it 


0.07 


1, 


,00 


2662 


1.52 ± 0.14 


1.31-12.00 


6 


0.30 ± 0.05 


1.31 ± 0.07 


1, 


,11 


153.14 ± 7.99 


-0.93 


± 


0.12 


0.67 


it 


0.07 


0.69 


it 


0.07 


1, 


,45 


2665 


1.99 ± 0.15 


1.42-8.00 


8 


1.22 ± 0.32 


1.42 ± 0.05 


0. 


,89 


106.03 ± 2.24 


-0.19 


± 


0.11 


0.75 


it 


0.12 


0.79 
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4. conclusions and discussions 

Pulses are the basic, central building blocks of GRB prompt emission, and it is essential 
to our understanding of GRB physics (Hakkila et al. 2008). In this paper we select a 
well-separated single pulses presented by Kocevski et al. (2003) and Norris et al. (1999) 
including bright and weak bursts to investigate the evolutionary slope of Ep during the pulse 
decay phase. We first fit the Ep of time-resolved spectra for our selected pulses with a single 
power-law function form and obtain the evolutionary slope to compare it with the theoretical 
predication of Paper I. Our results show that these observed evolutionary slopes are in good 
agreement with the predictions of Paper I. Contrary to Paper I, however, we find that within 
the limits of uncertainty the corresponding intrinsic spectra of most of bursts (about 25 for 
sample 1 and 29 for sample 2) may bear the intrinsic Comptonized or thermal synchrotron 
spectrum, whereas only a small number of bursts (about 16 for sample 1 and 12 for sample 
2) may be associated with the intrinsic Band function spectra. Paper I investigated 12 
pulses and found that there are eight consisting with intrinsic Band spectrum and only two 
corresponding to the intrinsic Comptonized or thermal synchrotron spectrum. We argue 
our analysis results are more rehable for the following reasons. Firstly, the quantity of our 
sample (56 for Band model and 59 for Compton model) much larger than that of Paper I. 
Secondly, we remove all possible improperly-measured Ep, i.e. the Ep with a < —2 and 
f3 > —2 for the BAND model, a < —2 for the COMP model are discarded, but Paper I does 
not seem to consider it. Lastly, the Ep with its error larger than 50 percent of itself are also 
excluded but Paper I includes the Ep with much larger error. Since the Paper I pointed out 
that the evolutionary curve of the decay portion of the pulse is dominated by the curvature 
effect our results also indicate that most of the FRED pulses we selected indeed result from 
this effect. 

There are also some outhers (about 15 for sample 1 and 18 for sample 2) included in our 
sample. Besides the curvature effect, there might be some other factors that can affect the 
value of S. One factor would be the variation of the rest-frame emission mechanism, which was 
revealed in Qin et al. (2005). For example, different rest-frame spectra or different speeds 
of the rest-frame spectral softening could lead to different values of the power-law index. 
Other factors pointed out by Paper I, such as the real intrinsic radiation mechanism may 
be more complicated as well as GRBs might be associated with more complicated situations 
rather than a simple fireball expanding isotropically with a constant Lorentz factor. There 
are some disparate Ep values that could be the result of source confusion (pulse overlapping) 
in the pulse decay phase (e.g. we can see from Figures 1 and 2). Maybe this is also one of 
the reasons giving arise to the outliers. 

We also study the possible correlations between the S and other parameters. A strong 
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the meanings of the parameters axe the same as table 1. 
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correlation between S and the Ep^^ax during the phase of pulses decay is identified. There 
might be some important link between S and the Ep^^nax- Also an important correlation 
between the observed Ep of pulse time-integrated spectra and S exists. In order to reveal 
if this correlation is intrinsic we attempt to use the pseudo-rcdshift calculated by Yonetoku 
et al. (2004) to check it. We also find there is correlation between the S and intrinsic peak 
energy Ep^i but for sample 2 the intrinsic correlation is looser than the observed one. One of 
reasons we consider is that the pseudo-redshifts may not well represent the real redshifts of 
our selected bursts. Another reason we suspect that the Band model may be more suitable 
for describing the GRB time-integrated spectra than Compton model. 

Ep is one of the fundamental characteristics of the prompt emission of GRB. Based on 
it several important empirical relations have been proposed. One of important relations is 
the so-called Amati relation (Amati et al. 2002; Amati 2003), i.e. the correlation between 
Ep in the GRB rest frame {E^^'^) and the isotropic radiated energy {E'^^°). The second 
correlation with much tighter correlation than Amati relation is between the E^'^ energy 
and the collimation-corrected energy (-E^), the so-called Ghirlanda relation (Ghirlanda et 
al. 2004). Similarly Liang & Zhang (2005) found a good correlation between EjJ''^, E^^°, 
and the achromatic break time in the afterglow light curve (tjet)- The third relationship is 
between E^^^ and the isotropic peak luminosity (I/*^"), the so-called the E^^^, L**° (Yonetoku) 
relation (Yonetoku et al. 2004). If these relations are vahd, they must be related to the 
fundamental physics of GRBs. Thus, Ep'^ energy provides us important knowledge about 
the characteristics of the prompt emission of GRBs. In order to calculate the bolometric 
fiuence which refiects the total radiated energy in the prompt emission the observed Ep is 
a crucial quantity. Therefore, the evolutionary slope might be a useful quantity to help 
us understand the physic of prompt emission since the evolutionary slope is related to the 
observed Ep as well as the intrinsic Ep. 

The observed Ep decay is clearly related to the pulse peak lag described by Hakkila et 
al. (2008) and Hakkila & Cumbee (2009). Both the Ep decay and the pulse peak lag describe 
pulse spectral softening from different viewpoints (the Ep decay reveals spectral changes with 
time, while pulse peak lag examines a self-similar temporal structure at different energies). 
There are several interesting and pertinent repercussions of this. 

Firstly, pulse evolution is not limited to the "pulse decay" phase. Although the pulse 
peak intensity is defined to be the time at which the decay begins, this is not the case: pulse 
peak lags indicate that the low-energy count rate is still increasing at the time of pulse peak 
intensity whereas the high-energy count rate has already started to decay. So, when does 
the pulse Ep decay begin if not at the time of pulse peak intensity? The evidence presented 
here from the plots in Figures 1 and 2 suggests that for most of pulses it probably starts 
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dropping at the pulse onset: the entire pulse likely represents a decay of Ep (this result has 
been mentioned before by various other authors, e.g. Norris et al. 1986). However, this 
observational result contradicts the theoretical conclusion of Lu et al. (2006), which states, 
"Now it is clear that, in the mechanism of the curvature effect, it is the decaying phase of 
a local pulse that contributes to the observed lags." Since Kocevski & Liang (2003) argued 
that the observed lag is the direct result of spectral evolution (the Ep decay through the four 
BATSE channels), something appears to be missing from the Lu et al. (2006) theoretical 
models, based on observations? We shall analyze it in detail in the future work. However, 
there are also some pulses whose spectra exhibit "tracking" (the observed spectral parameters 
follow the same pattern as the flux or count rate time profile) behaviors (e.g. Crider et al. 
1997; Peng et al. 2009). The trigger 5601 demonstrated in Figure 1 is such a "tracking" 
pulse. In this case, the Ep starts dropping at around the pulse peak intensity, while the Ep 
rises during the pulse rise phase. In fact the Hakkila et al. (2008) and Hakkila & Cumbee 
(2009) results have showed that spectral and intensity evolution simply scale with the pulse 
duration. Short, "tracking" pulses should thus exhibit the same type of Ep decay as the 
longer pulses, but vary on too short a timescale for this variation to have been recognized in 
many studies. 

Secondly, pulse evolution appears to be a defining characteristic of GRB pulses. The Ep 
decline correlates with a wide range of other pulse properties. Hakkila et al. (2008) showed 
that pulse duration is tightly correlated with pulse lag and anti-correlated with pulse peak 
luminosity, and it certainly appears that the same is true here (for example, the Ep decay 
correlates with the pulse decay time interval) . Similar to what is found in this work Hakkila 
& Cumbee (2009) also found that pulse peak lag anti-correlates with pulse peak fiux and 
correlates with pulse asymmetry, which reveals pulse shape seems related to the pulse decay 
rate. It seems to mean that pulse evolution is more than a characteristic of the pulse; it 
appears to DEFINE the pulse. 

Lastly, pulse evolution appears to drive the evolution observed in GRB bulk emission. 
Most of GRBs exhibit hard-to-soft evolution, this is due to the smeared measurement of 
overlapping pulses, and to the fact that hard, short, luminous pulses typically occur near 
the beginning of a burst. Our paper demonstrates that spectral evolution exists at the even 
finer level of pulse evolution. The fact that pulse properties are correlated (pulse peak lag 
dictates pulse luminosity, duration, hardness, and asymmetry) suggests that internal pulse 
evolution is responsible for pulse-to-pulse variations and ultimately for the burst evolution. 

In this study the pulses are limited to FRED pulses, which is based on the fact that 
the single pulses often been described by a fast rise and exponential decay as pointed out 
previously. Recently, Hakkila et al. (2008) and Hakkila & Cumbee (2009) suggested that 
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most pulses in long GRBs are really just variations of the same phenomenon. Therefore, our 
analysis results about evolutionary slope should apply to those non-FRED pulses. 
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